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The spectral theorem of many-body Green's function theory when there are zero 
eigenvalues of the matrix governing the equations of motion. 
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Hahn-Meitner-Institut Berlin, Glienicker Strafle 100, D-14109 Berlin, Germany 
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In using the spectral theorem of many-body Green's function theory in order to relate correlations 
to commutator Green's functions, it is necessary in the standard procedure to consider the anti- 
commutator Green's functions as well whenever the matrix governing the equations of motion for the 
commutator Green's functions has zero eigenvalues. We show that a singular- value decomposition of 
this matrix allows one to reformulate the problem in terms of a smaller set of Green's functions with 
an associated matrix having no zero eigenvalues, thus eliminating the need for the anti-commutator 
functions. The procedure is quite general and easy to apply. It is illustrated for the field-induced 
reorientation of the magnetization of a ferromagnetic Heisenberg monolayer and it is expected to 
work for more complicated cases as well. 

PACS numbers: 75.10.Jm, 75.70.Ak 



Many-body Green's function (GF) techniques make 
use of the spectral theorem to connect correlation func- 
tions to the Green's functions, which in turn are deter- 
mined from equations of motion that can be formulated 
conveniently in terms of a matrix equation. Should this 
matrix have zero eigenvalues, particular care must be 
taken, since they cannot be handled directly by the spec- 
tral theorem for the commutator GF, since it would re- 
quire a division by zero. The standard treatment of this 
case demands in addition to the commutator GF a knowl- 
edge of the anti-commutator GF (see references 0,0 and 
textbooks 3, 4J). An equivalent treatment in terms of the 
anti-commutator GF alone may formally appear to avoid 
this problem, but in practice one is faced with the same 
numerical difficulties. In this paper, we show that these 
difficulties can be overcome by using a singular-value de- 
composition (SVD) of the equations of motion matrix to 
eliminate the space connected with the zero eigenvalues 
(null-space), thus reducing the size of the problem and, at 
the same time, making the anti-commutator GF super- 
fluous. That is, one uses projection matrices furnished by 
the SVD to reformulate the spectral theorem in terms of 
a smaller number of commutator GF's whose associated 
matrix has no zero eigenvalues. We show that this pro- 
cedure is generally applicable and offers an alternative to 
the standard treatment; i.e. we show that a knowledge 
of the commutator GF's is sufficient. 

The appearance of zero eigenvalues in the matrix gov- 
erning the equations of motion is not an unusual event. 
We have met this problem in much of our previous work 
on the application of many-body GF's to the descrip- 
tion of thin ferromagnetic films @, IS II II II E3 We 
treated the Heisenberg ferromagnet including single-ion 
anisotropy and exchange anisotropy for several values of 
the spin S for single-layer and mult i- layer systems. A 
variety of decoupling approximations [111 Il2| to arrive 
at a closed set of equations of motion were considered. 
In most of these cases, the equations of motion matrix 



has at least one zero eigenvalue per layer and, if higher- 
order Green's functions Q are used to achieve an exact 
treatment of the single-ion anisotropy, there are two zero 
eigenvalues per layer. Calculation of many-layer systems 
therefore leads to matrices with a large number of zero 
eigenvalues. This, coupled with the fact that the matri- 
ces are in general nonsymmetric, leads in some cases to 
numerical difficulties which can best be circumvented by 
eliminating the null-space. The method for doing this 
may be considered as an alternative to the standard pro- 
cedure of handling zero eigenvalues. In the following, 
we first introduce our methodology and notation by re- 
viewing briefly the standard text-book treatment of zero 
eigenvalues, then we describe the new method which is 
based on the singular value decomposition of the ma- 
trix governing the equation of motion, and finally, by 
way of example, we apply the new method to the field- 
induced reorientation of the magnetization of a ferromag- 
netic Heisenberg monolayer. 

The standard procedure.— The equation of motion for a 
Green's function vector 



G, 



«4 



in energy space is (omitting the subscript lu) 

u((A;B)) v = ([A, B] v ) + (([A, H]_; B)) v , 



(1) 



(2) 



where (...) = Tr(...e~@ H )/Z is the thermodynamic ex- 
pectation value, H is the Hamiltonian and A and B are 
operators depending on the problem under consideration, 
and 77 = =pl refers to the commutator or anti-commutator 
Green's functions, respectively. 

Repeated application of equation (0 to the higher- 
order Green's functions appearing on the right hand side, 
(([A, H\-\ B)), results in an infinite hierarchy of equa- 
tions. In order to obtain a set of solvable equations, one 
has to terminate this hierarchy, usually by expressing the 
Green's functions of a particular order in terms of those of 



2 



lower order. In many applications, only the lowest-order 
functions are retained using a decoupling procedure of 
the form 



(([AH]-;B)) V ~ T((A;B)) V , 



(3) 



where T is the (in general non- symmetric) matrix ex- 
pressing the higher-order Green's functions in terms of 
linear combinations of lower-order ones. In compact ma- 
trix notation, the equation of motion is 



(ujI - r)G„ = A„, 



(4) 



where 1 is the unit matrix, and the inhomogeneity vector 
is A„ = ([A,B\ V ). 

It is now convenient to introduce the notation of the 
eigenvector method of reference |9( , since it is particularly 
suitable for the multi-dimensional problems in which 
many zero eigenvalues are likely to appear. One starts 
by diagonalizing the matrix T 



LTR = O, 



(5) 



where Jl is the diagonal matrix of N eigenvalues, 
lo t (t — 1, ...,N), Nq of which are zero and (N — Nq) 
are non-zero. The matrix R contains the right eigen- 
vectors as columns and its inverse L = R -1 contains 
the left eigenvectors as rows. L is constructed such that 
LR = RL = 1. Multiplying equation J3J from the left by 
L, inserting 1 = RL, and defining new vectors Q v = LG, ; 
and A n 



LA, ; one finds 



{uji - n)g v = A n . 



(6) 



Each of the components r of this Green's function vector 
has but a single pole 



(A v )t 



(7) 



This allows the spectral theorem, (see e.g. 0, HI) i to 
be applied to each component of the Greens's function 
vector separately: Define the correlation vector C — LC, 
where C = (BA) is the vector of correlations associated 
with G, ; . From the spectral theorem, the component r 
for the anti-commutator case (rj = +1) is 



(A) 



e /3u) T _|_ ^ 



(8) 



Use of the anti-commutator functions appears at first 
sight more straightforward (zero eigenvalues seem to 
make no problem), but leads in general to momen- 
tum dependent inhomogeneities, which are difficult to 
treat. However, it is usually easier to employ commuta- 
tor Green's functions, because in this case the inhomo- 
geneity vectors are momentum independent. 



For the commutator functions (r) = —1), the N — Nq 
components for ui T ^ 0, C 1 (the upper index 1 refers to 
the non-null space) , arc 



(-4 



VjT 



V 



(9) 



This equation cannot be used to define the No compo- 
nents C° corresponding to uj t — because of the zero in 
the denominator. Instead, one must enlist the help of the 
anti- commutator Green's function: 



(C°) T0 = lim - -w(^=+i) r . 



(10) 



The components of C , indexed by r (the upper index 
refers to the null space), can be simplified by using the 
relation between the commutator and anti-commutator 
inhomogeneities, Ai = A_i + 2C. This yields, for uj to = 
0, 



1 Hm uAA^ 

2 w^O UJ — LO Ta 

1 



= -(L (A_ 1+ 2C)) To = (L°C) n) , 



(11) 



where use has been made of the relation L°A_i = 0, 
which derives from the fact that the commutator Green's 
function is regular at the origin 0. 

At this point it is again convenient to introduce super- 
scripts and 1 to denote vectors corresponding to zero 
and non-zero eigenvalues, respectively. The right and 
left eigenvectors and the correlation vectors may then be 
partitioned as R = (R 1 R°) and 



L 1 
L° 



C = 



c 1 
c° 



(12) 



where the correlation vectors from equations © and i|ll|) 
are then C° = L°C and C 1 = £ l ~L l A, and S 1 is the 
(N - N ) x (N - N ) matrix with l/(e^- - 1) on the 
diagonal. 

Multiplying the correlation vector C from the left by 
R yields a compact matrix equation for the original cor- 
relation vector C: 



C = R^lM + R'VC, 



(13) 



This equation is in momentum space and the coupled 
system of integral equations obtained by Fourier trans- 
formation to real space has to be solved self-consistently. 
This was done in all our previous work |a, la la S lill • 

The singular value decomposition.— We now present 
an alternative method of treating the zero eigenvalues. 
Starting again from the equation of motion (0J , we write 
the r matrix as a singular value decomposition [lflj : 



uwv. 



(14) 
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The matrix W is a diagonal matrix whose elements are 
the singular values, which are > and U and V are or- 
thogonal matrices: UU = 1 and VV = 1, where V de- 
notes the transpose of V. U, V and W can be obtained 
very efficiently numerically 13]. However, one can also 
get the vectors V and U directly by diagonalizing IT or 
IT respectively. The singular values are just the positive 
square roots of the eigenvalues of these matrices, W 2 . 

vfrv = vvwuuwvv = w 2 , 
urfu = uuwvvwuu = w 2 . (15) 

Note that the T matrix is fully determined by the non- 
zero singular values and the corresponding eigenvectors: 



r = uwv = uwv. 



(16) 



u and v are N x (N — Nq) matrices obtained from U 
and V by omitting the columns corresponding to the zero 
singular values. The matrix w is the (N — No) x (N — Nq) 
diagonal matrix with the non-zero singular values on the 
diagonal. 

The crucial point is that the dimension of the equations 
of motion can be reduced by the number of zero singular 
values by applying the following transformations 



7 = vTv, 
g = vG, 
a = vA, 



(17) 



and realizing that T = Tvv and vv = 1. 

The (N — A^o)-dimensional equation of motion then 
becomes 



(ujI - 7)g = a. 



(18) 



Here again the eigenvector method can be used, where 
1 and r diagonalize the 7-matrix, l7r = or. Applying 
the spectral theorem to the equation of motion with the 
matrix 7 (which now has no zero eigenvalues) yields the 
equation for the correlations c = vC in momentum space 



c = rSha 



(19) 



where the (N—No) x (N—Nq) diagonal matrix E 1 has the 
same elements as before, \j(eP^ T — 1). In order to deter- 
mine the correlations in coordinate space, one has to per- 
form a Fourier transform and then solve self-consistently 
the system of integral equations 







dk^HvA- vC) 



(20) 



where A and C are the inhomogeneity and correlation 
vectors of the original problem, and the integration is 
over the first Brillouin zone. 



The advantages of the new procedure are that, after 
projecting onto the non-null space, it is no longer neces- 
sary to calculate the term in the spectral theorem which 
is connected with the anti-commutator Green's function, 
and the dimension of the final system of equations is re- 
duced by the number of zero eigenvalues of the original 
r matrix. 

Example: The Heisenberg monolayer.— As an example, 
we apply the new procedure to a ferromagnetic Heisen- 
berg monolayer with spin S = 1 , in which the magnetiza- 
tion can be reoriented by applying an external field. This 
is an example in which zero eigenvalues of the T matrix 
occur and is treated in references 0, 0] with the standard 
procedure not only for the S = 1 monolayer but also for 
the multilayer case and for S > 1. 

We restrict ourselves to a spin Hamiltonian consist- 
ing of an isotropic Heisenberg exchange interaction with 
strength Jki between nearest neighbour lattice sites, a 
second-order single-ion lattice anisotropy with strength 
Ki.h-, and an external magnetic field B = (B x , B v , B z ): 



U = - \Y.Wk S t+S z k Sf)-Y,K2,k{S z k ) 2 

(kl) k 

- E (l B ~ S k + \b + S-+B*s£). (21) 



Here the notation Sr: 



S% ± iS v k and B ± 



B x ± iBv 

is introduced, where k and I are lattice site indices and 
(kl) indicates summation over nearest neighbours only. 

Applying a Tyablikov (RPA) decoupling to the ex- 
change interaction term, and an Anderson-Callen decou- 
pling to the single-ion anisotropy term one obtains 0, 
equations of motion for the monolayer of the form 

(wl - r)G„ - A„ (22) 

with the Green's function and inhomogeneity vectors 



V 



A+.mn 

A™ n = \ A-' mn 

Az.mn 



G™"(k,c) = [ G-> mn (k,u>) 
G z v > mn (k,uj) 

(23) 

where G^™ n {u) = «Sf ; {S?) m (Sp n )) u , v and a = 
(+, — , z) takes care of all directions in space. For a field 
in the x-direction only, B = (B x ,0,0), the reorientation 
of the magnetization takes place in the irz-plane, and the 
r matrix for spin 5 = 1 (n > 1, m > 0, m + n < 2) is 









-H z 
\H X 



-H x 
H' 




(24) 



where 



H z = Z+(S z )J(q~J k ) , 
1 



H x 



K 2 2(S z )[l-—[S(S + l)-(S z S z }] 
B* + <S a! )J( g -7 k ) . 



(25) 
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Heisenberg monolayer: J=1 00, K 2 =1 , B x =0.4 
field-induced reorientation of the magnetization 




10 20 30 40 50 

temperature 

FIG. 1: Magnetizations (S z ) and (S x ) and the reorientation 
angle O for a spin S=l Heisenberg monolayer as function of 
the temperature. 



For a square lattice with lattice constant unity, 7k = 
2(cosfc x + cosfc y ), and q = 4 is the number of intra- 
layer nearest neighbours. The singular value decomposi- 
tion of equation (|14|l can be performed analytically for 
the monolayer case. The singular values of the diagonal 
matrix W are given as the positive square roots of the 
eigenvalues of the matrix IT or IT 



result, shown in figure 1, agrees perfectly with that ob- 
tained in references 0, Q obtained with the standard 
procedure, i.e. by solving equation 1)13(1 . In that work, 
spins S > 1 and multilayer systems were also considered. 
This can also be handled by the new algorithm but, in 
this case, the singular value decomposition and the ma- 
trix v have to be obtained from an appropriate numerical 
method El- 



Summary.— We have developed a new method for deal- 
ing with zero eigenvalues of the matrix governing the 
equations of motion for many-body Green's functions by 
applying a singular value decomposition to this matrix. 
This enables one to eliminate the null space, thus reduc- 
ing the number of integral equations which have to be 
solved self-consistently by the number of zero eigenval- 
ues of the equation of motion matrix. The calculation 
of the anti-commutator Green's function, necessary in 
the standard procedure using the spectral theorem, is 
superfluous in the new procedure, which is demonstrated 
for the example of the field induced reorientation of the 
magnetization of a spin S = 1 Heisenberg ferromagnetic 
monolayer. The method is expected to work for more 
complicated cases, where it has decided advantages over 
the usual procedure. 



W = 



ei 
e 2 




(26) 



where t x = ^H Z H Z + 2H X H X , e 2 = ^H Z H Z 
and one singular value zero. The matrices U and V are 
obtained from the eigenvectors of these matrices 



U 



-V2 
2 

2 





-H z 
V2e 2 
-H z 
V2e 2 

V2e 2 



Ml 

2(2 

Ml 

'-2 



(27) 



and 



V 



-H z 
V2ei 
-V2 

2 
H" 



V2ei 
2 



\/2H* 


Ml 

ei 



(28) 



The lower case rectangular matrices u and v introduced 
above are obtained by omitting the last column in U and 
the last row in V respectively. One sees explicitly that 
vv = 1 and that vv is a projection operator onto the 
non-null space because it has two eigenvalues equal to 1 
and one equal to 0, and (vv) 2 = vv. 

Having constructed the matrix v, we now solve equa- 
tion 1)20(1 for the spin 5 = 1 monolayer with the curve fol- 
lowing method, described in detail in reference Q. The 
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